Spectroscopy of low-frequency noise and its temperature dependence in a 

superconducting qubit 
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We report a direct measurement of the low-frequency noise spectrum in a superconducting flux 
qubit. Our method uses the noise sensitivity of a free-induction Ramsey interference experiment, 
comprising free evolution in the presence of noise for a fixed period of time followed by single-shot 
qubit-state measurement. Repeating this procedure enables Fourier-transform noise spectroscopy 
with access to frequencies up to the achievable repetition rate, a regime relevant to dephasing in 
ensemble-averaged time-domain measurements such as Ramsey interferometry. Rotating the qubit 's 
quantization axis allows us to measure two types of noise: effective flux noise and effective critical- 
current or charge noise. For both noise sources, we observe that the very same 1 / /-type power laws 
measured at considerably higher frequencies (0.2 — 20 MHz) are consistent with the noise in the 
0.01 — 100-Hz range measured here. We find no evidence of temperature dependence of the noises 
over 65 — 200 mK, and also no evidence of time-domain correlations between the two noises. These 
methods and results are pertinent to the dephasing of all superconducting qubits. 

PACS numbers: 03.67.Lx, 74.40.-n, 74.25. Sv, 85.25.Cp, 85.25.Dq 



I. INTRODUCTION 

A major remaining obstacle to implementing fault- 
tolerant quantum computation with superconducting 
qubits is the insufficient coherence time T2 compared 
to the gate-operation time. In the Bloch-Redfield pic- 
ture of two-level system dynamics, there are two mech- 
anisms that limit Ti\ energy relaxation T\ due to noise 
at the transition frequency Vq\, and dephasing T v due 
to low-frequency fluctuations of 2/01 • In the qubit's en- 
ergy eigenbasis, these contributions are due to the trans- 
verse and longitudinal components of the fluctuations, 
respectively. In cases when both relaxation and dephas- 
ing exhibit exponential decay laws, their inverse times 
add to a decoherence rate T-T 1 = {2T 1 y 1 + T' 1 . While 
Ti can now exceed 10 /is in several superconducting qubit 
modalities 1 - - — , further improvements are required to com- 
fortably exceed even the most lenient error correction 
thresholds. In general, energy relaxation is irreversible, 
such that further control-based improvements require 
resource-intensive multi-qubit quantum error-correction 
protocols. Dephasing, on the other hand, can be refo- 
cused by dynamical-decoupling techniques 2 -, with only a 
modest amount of resource overhead. The ultimate goal 
is to mitigate and eliminate the noise leading to both 
types of decoherence. To this end, a more detailed under- 
standing of the noise processes - such as magnetic-flux, 
critical-current, and charge fluctuations - would expedite 



materials science, device engineering, and the develop- 
ment of coherent-control methods. 

Effective surface spins have recently been identified 
as one dominant source of low-frequency magnetic-flux 
noised, detrimental to several types of superconducting 
qubits; however, open questions remain regarding the na- 
ture of these spins. Their noise is known to be due to local 
fluctuators 6 - - — and the spectrum exhibits a 1//" power- 
law dependence from hertz to tens of megahertaifir— . Its 
dependence on the device geometr y 15 ! 16 merits further 
study. 

Similarly, for the flux qubit 17,18 , the noise in the tunnel 
coupling A between the persistent-current states shows a 
l//-type spectrum from hertz to hundreds of kilohcrtz 2 . 
This noise may originate in the critical-current fluctua- 
tions of the Josephson junctions 1 ^— , and/or fluctuating 
offset charges, due to, e.g., charge traps located in the 
oxides of the junction, metal-insulator interfaces, or sur- 
faces. Charge noise can lead to dephasing even in the 
flux qubit, even though the junctions have a relatively 
high ratio of Josephson-tunneling to Coulomb-charging 
energies (Ej/Ec ~ 50 in our device). 

In this paper, we introduce a measurement technique 
for low-frequency noise. A distinguishing feature of our 
technique is that it enables the measurement of noise 
spectra up to frequencies limited only by the achiev- 
able measurement repetition rate. This is important, 
because noise measured in this manner resides (at least 
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in part) within the relevant measurement bandwidth of 
time-domain experiments, e.g., Ramsey interferometry, 
that use the standard ensemble-averaging (that is, the av- 
eraging of multiple trials acquired at the same repetition 
rate) to estimate the qubit-state occupation probability. 

We report on a direct characterization of the 1 / /-noise- 
power spectral densities (PSD) S(f) in an aluminum su- 
perconducting flux qubit^. We distinguish between the 
two noises Se, which is effective flux noise, and SA, which 
can be parameterized as effective critical-current noise 
or effective charge noise. Interestingly, we find that the 
same 1// Q power laws, measured at much higher frequen- 
cies^, extend down to the 10~ 2 — 10 +2 -Hz range nearly 
unchanged. Over the temperature range 65 — 200 mK, 
both noises, Se and 8 A, are independent of temperature, 
and any Se — <5A-noise correlations are very small or non- 
existent. 



II. EXPERIMENTAL METHODS, ANALYSIS, 
AND RESULTS 

The flux qubit's two- level Hamiltonian is H = 
— (h/2)(ea x + Aa z ). Here e = 2I p $b/h is the energy 
detuning between the diabatic states of classical circu- 
lating current I p = 0.18 ^A, and e is adjusted by the 
external magnetic flux $ via <£>b = $ — ( 3 ) o/2 ($o = h/2e 
is the superconducting flux quantum); see the schematic 
in Fig. QJa-b). A = 5.4 GHz is the tunnel coupling 
that hybridizes the persistent-current states, and is es- 
tablished during fabrication by the Ej/Eq ratio. We 
write each of the parameters A = e, A as the sum 
of its nominal value and a time-dependent fluctuation, 
X(t) = A(°> +SX(t). We distinguish between the effects of 
Se and SA fluctuations by rotating the qubit's quantiza- 
tion axis (eigenbasis) , thereby altering the sensitivity of 
the energy-level splitting v m = \J e 1 + A 2 to fluctuations, 
D\ = dvoi/d\ — A/Voi; see Fig. QJb). The dominant 
contributor to longitudinal fluctuations in the qubit's en- 
ergy eigenbasis is SA noise at e = (as the second-order 
contribution from Se noise is negligible 2 -), and Se- (flux) 
noise for |e| > 0.1 GHz. 

We use a hysteretic SQUID to read out the qubit's 
state. During qubit manipulation, we conveniently set 
the SQUID current lb to zero. (This is close to its op- 
timal value of 1^ — 50 nA in this device, at which SI\, 
fluctuations are decoupled from the qubit 2 ^.) 

To directly probe the fluctuations, we repeatedly let 
the qubit undergo Ramsey free induction; see Fig. HJc). 
Instead of scanning the pulse separation r, we fix it 
at a value To, chosen to maximize the qubit state's 
sensitivity to noise. We also fix the nominal detun- 
ing z/°) of the applied microwave frequency from the 
qubit's frequency to = vq\ — — (a free- 

evolution 7r/2 rotation in the X—Y plane); see Fig. [TJd- 
f). The SX(t) fluctuations translate into frequency fluc- 
tuations, v(t) = + Sv(t), which in turn translate 
into fluctuations of the SQUID's switching probability, 




Pulse separation t (us) 

FIG. 1: (color online) (a) Schematic of the qubit device with 
biasing and read-out circuitry, (b) Qubit energy-level dia- 
gram near <!>b = $o/2 (e = 0). (c) Ramsey interference 
fringe measured at e = with 5, 000 averages per point, (d) 
Double-7r/2 pulse (Ramsey) sequence, (e) Bloch sphere rep- 
resentation of dephasing and our measurement scheme: the 
initial 7r/2 pulse takes the Bloch vector from the ground state 
(along Z) to the equator (along Y); during the free-induction 
time t = ro, the nominal detuning is = l/4ro, but differs 
slightly at each trial (quasi-static noise), and the Bloch vec- 
tors therefore acquire different phases; the final 7r/2 pulse ro- 
tates the Bloch vectors out of the equatorial plane so that they 
can be distinguished by the SQUID read out. (f) Calibration 
by scanning i/ MW with ro = 0.3 fis, e = 0. The red, solid line 
is a sinusoidal fit used to determine the working point for 
maximal sensitivity (dashed line at z/ ' = —0.8 MHz), (g) 
Individual SQUID-switching events measured on an oscillo- 
scope. The dashed, horizontal line is a software threshold 
detector. 



Psw{v, t ) = p — a (T) cos(27r v{t) t ), where a (T) is the 
temperature-dependent read-out visibility (2 ciq — 79 % 
at the refrigerator's T = 12-mK base temperature) and 
Po is the switching probability for the qubit's 50-% su- 
perposition state. We linearize about a working point at 
Psw = Po and obtain dp sw /dv = 2tt eto(T) tq. Due to this 
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Frequency f (Hz) Frequency f (Hz) 

FIG. 2: (color online) Bilateral noise PSDs. (a) Standard-method PSD of the ensemble-averaged time series of switching 
probabilities p S w (blue data, on left) and of the binary time series of single-shot measurements {z n } (gray, on right), measured 
at base temperature. The red points (thin line) are averages of hundreds of traces. (S E measured at e — 450 MHz.) (b) 
Cross-PSDs of interleaved time series, S\(f) (Eq. 0, vs. device temperature (S e measured at e — 640 MHz). The data were 
smoothed by a sliding average with a triangular weight function of width A/ = //4. The dashed, diagonal lines are the 1/ / 
noises S\(f) = (2iv) 2 K\Ai,j c /\f\, derived in Ref. y; the sensitivities are ka,i c = dA/di c ~ 3.2 GHz (i c is normalized critical 
current) and k s = de/d<& = 1.1 GHz/m$o, and the noise strengths are Ai a = (4.0 x 10~ 6 ) 2 and = (1.7/i$o) 2 . The solid, 
diagonal line is S e (f) = (27t) 2 k 2 A^/|/| ' 9 , with ^44 = (0.8 /i<J>o) 2 - See Ref. 2 and its supplementary information for details. 
The horizontal, green, dashed lines are the sampling- noise levels Sn(ro, T) at low temperature; the triangles represent an upper 
cut-off frequency f c for sufficient averaging, above which the data are not dependable. 



transfer function, a correction factor arises in the calcula- 
tion of ensemble-averaged quantities, assuming Gaussian 
statistics (see section B.2 of the appendix). This fac- 
tor is a(ro,T)/ao{T), where cl(to,T) is the amplitude of 
the fringe at pulse separation tq- The conversion factor 
from the noise i5A to the switching probability p sw then 
becomes ijx (r ,T) = 2tt a (to, T) t D\. 

The single-shot read-out of the qubit, with a repeti- 
tion time At, results in a binary time series {z n }; see 
Fig- Dig)- (I n thfs experiment, At = 2 ms to allow for the 
read-out induced quasiparticles to relax between trials so 
that they contribute negligibly to heating.) Each element 
z n represents the result of a Bernoulli trial with expec- 
tation value p sw . The standard method to determine the 
noise-PS D 13 ' 15 i 24 is to ensemble average the switching 
events acquired during a gate time t acq = N(jAt, with 
typically Nq ps 1,000, to determine the average switch- 
ing probability p sw (binomially distributed), and then 
take the Fourier transform of the time series of switch- 
ing probabilities p sw ; see Fig. (Ha). This approach esti- 
mates the noise spectrum at frequencies up to l/2i acq . 
In contrast, in our approach, we calculate the bilateral 
PSD from the recorded series of single-shot measure- 
ments (Bernoulli trials), S(f k ) = \Z k \ 2 /(N/At), where 
fk = k/NAt, k = 0, . . . , N/2, and {Z k } is the discrete 
Fourier transform^ of {z n }, with n typically ranging 
from 1 to N = 5 x 10 4 . This method increases the upper 
cut-off frequency from l/2NoAt to 1/2 At, which may 



approach 1/to and is limited only by the achievable rep- 
etition rate; see Fig. [2ja-b), section A of the appendix, 
and Fig. 

Both the PSDs originating from single-shot 
and from ensemble-averaged measurements 
in Fig. [Ha) exhibit statistical-sampling noise 
S n (r ,T) = (27r) 2 cr s 2 At/ril(T ,T), where the vari- 
ance is <7 2 = p sw (l — Psw)- We can eliminate this 
white background noise by calculating the bilateral 
cross-PSD of the two interleaved (single-shot) time series 
z' = z 2n -i and z" = z 2n , 



S CIOSS (f k ) = 



Z' k (Z' k T 
N/At ' 



fk = 



NAt' 



(1) 



where now k — {),... , N/4. Dividing Eq. ((T|) by the con- 
version factor, we obtain the spectral density of the fluc- 
tuation A, 



S x (f k ) = (2nyS c ™ s %f k )/r 1 i(T ,T). 



(2) 



We typically average the spectra of 500 time series to im- 
prove statistical accuracy without compromising band- 
width, and recalibrate the working point periodically 
(hourly). Note that both the 1/f noise and the sam- 
pling noise dominate all other background noise at the 
temperatures considered. 

The Se and SA noise-PSDs are plotted in Fig. E^b) 
for several temperatures. There is striking agreement 
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with the 1//" power laws inferred in Ref. |U, measured at 
considerably higher frequencies (0.02 — 20 MHz). Noise 
that is strictly 1// Q=1 over the frequency range rel- 
evant to free-induction (here 10 _1 ~ 10 6 Hz) gives a 
Gaussian decay function of the temporal Ramsey os- 
cillations. Assuming that our SA noise satisfies these 
criteria, we use the approach of Ref. [25| to calculate 
the inhomogeneously broadened decay-time constant T*. 
With noise sensitivity ka,i c and strength A; c as de- 
fined in the legend of Fig. I\ = 12 /is, and typi- 



cal parameters r 
T 



1 MS and i acq = 



- (2ir KA . lc D A )- L A ic 1/2 (ln[i acq /2T ])- 1/z = 3.2 (is, 
in very good agreement with the observed T* in Fig. (He). 

We now turn to possible Se — <5A-noise correlations in 
the time domain, to check that the two spectra in Fig. [2] 
are not due to one and the same mechanism. Figure [3] 
shows how we repeatedly measured the switching proba- 
bility at alternating flux biases, p sw (±e^+Se, A^+SA), 
with ±eW chosen such that the effects of the two noises 
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on Sv were similar in magnitude, i.e. D e (e 



(0) 



)\Se 



Da ( £< ^) We set the pulse separation tq and nomi- 
nal frequency detuning With the energy-level split- 
ting u(e, A) = ( e (°), A(°)) + Sv (e(°), Se, SA), we use 
the decay function of the Ramsey fringe to infer the noise 
correlations from the measured p sw , 
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p sw (e, A) = po - a exp (-to/2T x ) 



(3) 



x exp 



(- [r /r;(e)] 2 )cos (2nv(e,A)T ), 



where [l/T*(e)] 2 = [1/T*(0)} 2 + KA^D 2 and K is a con- 
stant that we have determined independently, along with 
the other parameters in the equation. At the bias points 
e = ±e(°\ Se fluctuations induce negatively correlated Sv 
fluctuations (and consequently p sw fluctuations), whereas 
SA fluctuations induce positively correlated Sv fluctua- 
tions. At each time step, the measurement of p sw for 
e = ±£o yields a system of two non-linear equations in 
the two unknowns Se and <5A. We solve this system nu- 
merically: Fig.[3ta) shows the rawp sw data and extracted 
Se and SA vs. time. We then calculate the cross-PSD 
SsA(f) as the Fourier transform of the cross-correlation 
function, and obtain an upper bound on the magnitude 
of the correlation function, as shown in Fig.^b), 



|7 £ a(/)| 



\SeA(f)\ 



[S e (f)S A (f)]^ 



< 0.2. 



(4) 



Finally, we measured the temperature dependencies of 
the two types of noise. Figure 2] shows the integrated 
noise powers I1a vs. temperature T in the 65 — 200-mK 
range, where our read-out visibility is sufficient. We 
observe in essence temperature independence for both 
noises. For Se (flux) noise, this is consistent with previ- 
ous observations in SQUIDs 5 .^; we discuss the 5 A noise 
below. 



FIG. 3: (color online) Se — <5A-noise correlations measured at 
base temperature, (a) Top panel: Time series of p sw with 
fixed r = 0.3 us, nominal v m = -0.8 MHz, and 1,000 av- 
erages per measured point. The flux-bias polarity (±£^) 
was alternated between positive (blue) and negative (red) 
[±$b = ±0.058 m«l>o, see inset in (b)], with a 2-s repetition 
period, (a) Middle and bottom panels: Differential (Se) 
and common-mode (SA) noise inferred from the data in the 
top panel, (b) Magnitude of the correlation function, (Eq.p4|; 
the smaller the magnitude, the more indeterminate the phase, 
and therefore we cannot discriminate between possible corre- 
lation and anticorrelation. Inset: 30-min measurement of 
Ramsey fringes (frequency detuning referred to e = 0, where 



voi = A 



White circles indicate the bias points used 



in the top panel of (a). 



III. DISCUSSION 

In order to analyze the i5A noise, we parameterize it 
as an effective, normalized critical-current noise, i c = 
SI C /I C , with I c = OA/iA, in a Josesphson junction with 
area A = (0.2 /im) 2 . Van Harlingen et alJ^ found a 
"canonical" value for the 1// <5/ c -noise power at 1 Hz 
and 4.2 K: A c £ n w 144 (pA) 2 (7 c /^A) 2 /(„4/pm 2 ) in sev- 
eral SQUIDs and qubits of various sizes, made of dif- 
ferent materials. The authors hypothesized a quadratic 
temperature dependence, consistent with certain plausi- 
ble models for the noise sources below 100 mK, while not- 
ing that other models suggest a linear dependence. The 
bilateral normalized noise-PSD then becomes <S' J can (/) = 
A c £ n I~ 2 (T/4.2K) 2 /|/|, which, for T = GSmK^is con- 
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FIG. 4: (color online) Temperature dependencies of the to- 
tal noise powers n A (T) = (27r) 2 K i / dfSx(f;T) in the fre- 
quency intervals F e = 0.02 - 50 Hz and F A = 0.02 - 2 Hz, 
cf. Fig. HJh). It is possible to measure the Se noise up to 
somewhat higher temperatures and frequencies than the 5 A 
noise. Note that the H\ values depend on the integration 
limits although the choice of Fa does not make any signif- 
icant difference in the trends. (The double data points for 
the 165- and 180-mK 8A noise were measured with differ- 
ent pulse spacings to.) We estimated the device temperature 
from switching-current measurements on the SQUID, which 
suggest saturation at dilution-refrigerator temperatures be- 
low about 65 mK. The error bars are derived only from the 
fit error of the read-out visibility a(ro,T), included in the 
p s „-to-<5A conversion factor r)\(To,T), see Fig. [Tf f). 



In conclusion, our spectroscopy of both Se noise (flux 
noise) and SA noise (effective critical-current or charge 
noise), facilitated by single-shot measurements and thor- 
ough data analysis, shows that the very same l// a de- 
pendencies, measured at substantially higher frequencies, 
extend down to millihertz frequencies. This apparently 
indicates that the same noise mechanisms are active and 
dominant over some ten orders of magnitude or more for 
Se noise and at least eight orders of magnitude for SA 
noise. The Se noise may extend, with roughly constant 
slope (on a logarithmic scale), up to the qubit's transi- 
tion frequency at several gigahertz^: there, this noise is 
nearly transverse to the flux qubit's energy eigenbasis, 
and would therefore also contribute to energy relaxation. 
The small, if not negligible, Se — 5A-noise correlations 
(over 7 x 1CP 4 — 2 x 10 _1 Hz) show that the noises are 
due to distinct underlying mechanisms. Moreover, both 
noises are temperature independent in the 65 — 200-mK 
range, which suggests that the microscopic mechanisms 
are dominated by even lower energy scales than that. 
This is useful information for the development of noise 
models. It also calls for further studies of the repro- 
ducibility of the device properties, and, in particular, of 
the SA noise, as it limits the coherence time in supercon- 
ducting flux and transmon qubits. 
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siderably lower (almost 20 times) than our measured 
value. On the other hand, Eroms et ali22 measured re- 
sistance fluctuations in aluminum tunnel junctions: they 
found about 100 times lower noise power at 4.2 K, a lin- 
ear temperature dependence, and saturation below 0.8 K, 
i.e., S ic (f) = (1/100) x A^ n J- 2 (T/4.2K)/|/|. With 
T = 0.8 K, this gives a value about 2.5 times lower than 
what we observe. We also note that recently, contrary to 
these findings, Paik et al. 3 reported no evidence for 1/f 
i c noise in a Josephson junction. 

An alternative source of SA noise is the fluctuat- 
ing offset charges, SQ, known to exhibit 1/f noise^— ; 
these charges effectively supply a gate voltage to each 
island. The charge-noise power typically observed in 
single-electron tunneling (SET) devices is proportional 
to temperature 2 ^ (although quadratic dependencies have 
also been observed 2 ^) and saturates below about 200 mK, 
due to self heating of the SET, at a "canonical" 
value Aq of about (1 ~ 10 me) 2 at 1Hz. We esti- 
mate our qubit's maximum sensitivity to charge fluc- 
tuations, ka,q = dA/dQ, to be in the range 0.1 ~ 
lMHz/e. We can then parameterize the SA noise as 
charge noise and estimate the dephasing time T* = 
(2^ KA , Q Z? A )-M Q 1/2 (ln[t acq /2r ])- 1 /2 „ 4 ~ 400 /js. 
The lower end of this range is not far from our observed 
value. Moreover, the tunneling of charged quasiparticlcs 
between the small islands constituting our device may 
displace offset charges and contribute to dephasing at 
e = 0. 



A. Spectral density and the statistical noise floor 

Here we describe how we calculate the noise-power 
spectral density (PSD) from the noisy time series, and 
eliminate the statistical white-noise floor due to sam- 
pling. 



1. PSD 

The fluctuations of our qubit's transition frequency 
constitute a zero-mean, wide-sense stationary process 
Sv(t); at our chosen working point, 8v(t) = v[t) — 1/4-rn. 
We seek its bilateral noise-PSD (in units of rad/s, i.e., 
we use the angular-frequency correlator), 



S(f) = (2ir) 2 lim 



<M/)I 2 > 



T 



(5) 



Our measurements' raw data, however, consists of a bi- 
nary time series {z n } with elements of expectation value 
Vn = (z n ) = Pa + o-o sin(27r Sv(nAt) To), where At is the 
time step. The statistical properties of {z n } represent 
those of the underlying process 5v{t), up to a conversion 
factor and a correction factor (explained in the next sec- 
tion). We can therefore take the discrete Fourier trans- 
form {Z k } = F[{z„}], identify Z(f k ) = Z k x At for f k = 
k/NAt, and compute the discrete, bilateral noise-PSD 
over the frequency range from 1/f tot = 1/NAt — lOmHz 
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to l/2t acq = 1/2 At = 250 Hz, 



S k = = (27T) 1 



2 iVAt ' 



5 fe#0 = (2tt) ; 



^ fc 2 (At) 5 



NAt 



(6) 



for fc = Q,...,(iV-l)/2 with / fc 



NAt 



We then take the statistical average of M different PSDs 
obtained from different time series measured in succes- 
sion, 



1 M ) 

m— 1 



(7) 



pulse height, and the stochastic variable ti is the pulse- 
arrival time. The Fourier transform of w(t) is W(f) = 

G(/)Ef=i^exp(-i27rA), where G(f) = F[g(t)]. Car- 
son's theorem is then 



{fk) = (2nf «l/At)(6 2 }|G(/ fe )| 2 + (b) 2 6(f k )) 



5 Carson 
n 

(9) 

In our case, the pulse height bi is binary and the pulse- 
arrival rate is fixed at 1/At; we can therefore write (with 
Kronecker's delta) G = J*[<S m ,o] = 1 x At. We just have 
to replace the mean-square (b 2 ) by the variance a 2 and 
set S(fk) = 5 k ,oAt to obtain Eq. ([5]). 

Parenthetically, one can also derive Eq. ([5]) by using the 
Wiener-Khintchine theorem. The autocorrelation func- 
tion is 



and finally smooth the result with a sliding average in 
the frequency domain. 

With this method (Eqs. EH?]), each element z n is the 
result of a single-shot measurement; the sampling time 
step t acq is the same as the pulse-sequence repetition time 
At. This sets it apart from the standard approach of first 
taking the ensemble average of typically Nq = 1,000 
samples in the time domain, before calculating the PSD 
of the resulting N/Nq sampled points. The acquisition 



N-l 



time is then t a 



iVcAt, and the upper cut-off fre- 



quency becomes only l/2t acq » 0.25 Hz. 



2. White-noise floor 

The PSD of the single-shot time sequence suffers from 
statistical sampling noise because each time step n consti- 
tutes a Bernoulli trial (b n ): the read-out SQUID switches 
(b n — 1) with probability p and does not switch (b n = 0) 
with probability 1 — p. This statistical noise has a white 
spectrum; it dominates possible white background noise 
from other sources, and dominates also the 1// noise at 
high frequencies. To estimate it, we can treat the stochas- 
tic variable b n as independent and identically distributed 
(i.i.d.) with ensemble-averaged mean (b) = p and vari- 
ance a 2 = ((Ab) 2 ) = (b 2 ) - (b) 2 = p(l-p). Sampling 
at a fixed rate 1/At, the white-noise floor of the bilateral 
PSD becomes 



S n (fk) = (2tt) 2 (a 2 + (b) 2 5 kfi ) At. 



(8) 



Here we use Kronecker's delta S k ,o in the discrete PSD. 

The same expression is valid for the PSD of the 
ensemble-averaged time series, the constituent elements 
of which have a binomially distributed switching proba- 
bility (c for "counts") averaged over a gate time t acq = 
NcAt : we obtain Eq. (|5|) after substituting (c) = p and 
ct 2 = p (1 — p) /Nq for (b) and a 2 , respectively. 

Equation ([8]) is, in fact, a modification of Carson's 
theorem, which is valid for temporally random pulse 
arrivals. There, one considers a random pulse train 
w(t) — YaLi big{t — ti), in which g(t) is the pulse 
envelope, the stochastic variable 6/ is the (continuous) 



R { ™ ] = (1/N) x Kbn- m = (b) 2 + o 2 8 m , , (10) 



and J-"[l] = Sk.o X At, so that 



s n (fk) = (2tt) 2 t kr' 



= Eq. 



(11) 



We find that the PSD resulting from a simulation of 
Bernoulli- and binomially distributed noise agrees well 
with the measured data and with Eq. ([5]): we therefore 
conclude that our experimental noise floor is due to the 
statistical sampling. 

If the data consisted of a train of pulses of finite length 
in time, the PSD would have a roll-off near the Nyquist 
frequency l/2At. For example, the Fourier transform 
of a boxcar (square) pulse of length At is the func- 
tion At sine (nf At). In our case, after conversion of the 
SQUID's response (the presence or absence of a voltage 
pulse) to binary form, our data can be seen as represented 
by a train of delta-functions, and their Fourier transform 
is frequency independent, i.e., our white noise floor has 
no roll-off. 



3. Cross-PSD: white-noise elimination 

In order to eliminate the white- noise floor, at the ex- 
pense of a halved Nyquist frequency, we calculate the 
discrete cross-PSD of interleaved time series, i.e., by set- 
ting z' n = z 2n -i and z" = z 2n (with n= I,..., N/2) and 
computing the cross spectrum of z' n with z". (We again 
assume that the stochastic switching process is uncorre- 
cted from sample to sample, at frequency 1/At.) The 
resulting PSD is 



ocross (c\ \A 

i> k = - W 5 N/2At 



s: 



fe#0 



= (2^ 



N/2At ' 
(12) 

where k = 0,...,N/4 and f k = k/NAt. These ex- 
pressions reproduce the noise spectrum, with the use of 
correction factors, as explained in the following section 
(Eqs. E 
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Compared to the previous section, we have eliminated 
the white noise by circumventing the zero-delay autocor- 
relation term in Eq. (fT0|) . and are only left with the 
delta-function component, 



Sr ss (fk) = (2ir) 2 (b) 2 5 kfi At. 



(13) 



In the same way as in Eqs. (fTUHTTj) , this can be de- 
rived by applying the Wiener-Khintchine theorem to 
the cross correlation function, which this time simply 
gives R^yyi — (b) 2 (the subscript b'b" indicates two in- 
terleaved sub-series obtained from the original series b n ), 
and therefore 



Sr ss (fk) = (2n) 2 T R 



Eq. (J 



(14) 



When calculating the PSD, we take the statistical aver- 
age of S"° ss , keeping the averaging coherent throughout 
(i.e. retaining S" oss as a complex quantity), and, just as 
for Sk, smoothen it with a sliding average before plotting 
its magnitude | (<^ ross ) stat >f q | . 

The result (|T2"j) is equivalent to the explicit subtraction 
of the incoherent noise S n from St (Eqs. [SI 0), with the 
advantage, however, of drastically reduced uncertainty, 
in particular at high frequencies where the (1 //-noise) 
signal is much smaller than the white noise. This method 
is appropriate for the analysis of, e.g., l//-type noise. 
However, it is not applicable in a predominantly white- 
noise environment: then, the noise under study would be 
eliminated along with the statistical white noise. 



B. Correction factors: Quasi-static noise and the 
non-linear transfer function 

In this section, we treat the effects on the PSD caused 
by quasi-static noise, and by the sine nonlinearity in the 
conversion from the measured switching events to the 
variations of the qubit's transition frequency. 



Decay of the Ramsey fringe - quasi-static noise 



The higher integration limit is here the inverse of the 
free-induction time, 1/r w 0.1 — 100 MHz; fluctuations 
at even higher frequencies are effectively canceling out. 
The lower limit is given by the total acquisition time i acq 
used to infer the qubit's population at each fixed free- 
induction time span r. Typically averaging over -/V avg = 
5,000 measurements with a repetition time At = 2 ms, 
we obtain i acq = N avg At = 10 s. (If instead the measure- 
ments were done in the opposite order, stepping over r 
in the inner loop, with N pts ~ 100 steps, and averaging 
over N avg in the outer loop, the total acquisition time 
would be N pts N avg At = 1,000 s, and the lower cut-off 
frequency would be correspondingly lower.) 

Ensemble averaging over all realizations of 5<p(r), and 
assuming Gaussian fluctuations resulting from numerous 
fluctuators, we obtain the dephasing envelope 

h(r) = (exp(i^(T))) = exp(-((M 2 >/2) = (16) 



exp 



/tog 



d/S(/)sincV/7 



« exp(-CT 2 r 2 /2) , 

where the sinc-squared function is due to the square time 
window of the Ramsey pulse sequence, and we can ap- 
proximate it by unity for / < l/2r. 

As an illustration, we now evaluate h(r) for the two 
cases of 1// noise and white noise. For 1// noise, S(f) = 
A/f, Eq. jUJ becomes a 2 = 2 A log [(l/2r)/ (l/t ac q)]- 
The weak, logarithmic sensitivity to the cut-off frequen- 
cies effectively allows us to treat it as a time-independent 
constant, a 2 ~ 2 AC, giving Gaussian decay pi/f(r) = 
cxp( — o~ 2 t 2 /2). For white noise, S(f) = S w , on the other 
hand, the integral is linearly sensitive to the upper cut- 
off frequency, so that a 2 — S w /t, yielding an exponential 
decay p w (t) = exp(— S w t/2). Here the exponent is pro- 
portional to time; we can therefore identify 1/T* = S w /2 
as the dephasing rate. 



Noise in the effective longitudinal field coupled to the 
qubit results in decoherence of the quantum superposi- 
tion. We denote a fluctuation as "quasi-static" or "inco- 
herent" noise, when it can be considered as static during 
each free-induction period, but varying over the longer 
time span between experimental realizations. Dephasing 
results from such uncorrelated fluctuations of the Lar- 
mor frequency vox, and therefore of the accrued phase of 
the superposition state, <p(r) = 2tt j dt i/oi(t). It leads 
to decay of the Ramsey free-induction signal, as each 
measured point is the incoherent average of many exper- 
imental realizations. We describe this fluctuation by a 
standard deviation, 

a 2 = 2 f 1/T dfS(f). (15) 



2. Repeated fixed-time free-induction 

The previous section described how quasi-static noise 
determines the dephasing of the Ramsey-fringe. Now we 
turn to its effect on Ramsey interference with a fixed 
free-induction time To, repeated numerous times. 

With our single-shot measurements, each element of 
the binary time series {z n } is a Bernoulli random vari- 
able z n with expectation value given by the switching 
probability p sw , which we now denote as 

Vn =Po + a smx n . (17) 

This function has a non-linear dependence on x n = 
27r<5!/„To, the phase accrued during To, where 8v n is the 
average fluctuation of the transition frequency at time 
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Resolved by Range relevant for 

1 ensemble avg. 1 time-domain meas. 
M »-H 

Resolved by 




tot acq (avg) acq (sing) 



FIG. 5: Sketch of the PSD, indicating the frequency inter- 
vals resolved by the ensemble-averaging and single-shot sam- 
pling methods. Also indicated are the Gaussian, quasi-static 
noise and the variances a\ 2 (Eq. 1 1 5p . Here t to t is the total 
length of the time trace (can be several minutes to hours); 
iacq(avg) = A^At = 1 ~ 10 s is the acquisition time per mea- 
sured point in time-domain experiments such as Ramsey and 
spin-echo decay; t a cq(sin g ) — At = 2 ms is the repetition time 
(acquisition time of the single-shot samples); and to as 1 fis is 
the pulse spacing. 



step n. This phase x n , in turn, has noise contributions 
from two distinct frequency intervals, "1" and "2." 

We denote as "interval 1" the frequencies which we 
can resolve by taking the Fourier transform of the series 
{z n }, of total length NAt and step size At, i.e. from 
l/t tot = 1/NAtxs lCT 2 Hz to l/2t acq = 1/2 At « 250 Hz 
(or with the interleaving method up to l/4At ps 125 Hz). 
The noise within this interval has zero mean and variance 
a\ (Eq.HSJ. 

In addition, there is a contribution from the quasi- 
static noise in "interval 2," which is the range from 
l/2i a cq to 1/t ; see Fig. This noise cannot be re- 
solved, but acts in aggregate and leads to dephasing, e.g. 
in a Ramsey- fringe experiment. It has zero mean and 
variance o\ (Eq. [15]) . 

Noise at even higher frequencies than 1/ro averages out 
during free induction. 

At each time step n, the element x n is subject to 
noise contributions from both intervals, and their vari- 
ances add up to (72 = of + <j\. We write x n = u n + v n , 
where u and v refer to the noise originating in intervals 
1 and 2, respectively. Here u n has correlations between 
the different time steps n due to the memory effect of the 
1// noise; on the other hand, v n is incoherent and can 
be taken as a Gaussian i.i.d. random variable. 

While it is impossible to unequivocally infer x n from 
the measured z n at each instance n, we can infer statisti- 
cal properties of {x n }, such as its correlations and spec- 
tral density, up to the frequency l/2i acq = 1/2 At, which 
can approach 1/tq. This is advantageous compared to 
the ensemble-averaging method, which has a longer ac- 
quisition time i acq = N&At. 

We can write the m ^ n autocovariance function for 



Az n — z n (z n ) 8-S 

(Az m Az n ) = (Ay m Ay n ) = ao(sina; m sina;„) ~ (18) 

— Q,Q^X m X 7 /) — Q,Q^U rn Un) • 

The first equality holds because the Bernoulli trials are 
independent, and the last equality is the consequence of 
v n being i.i.d., which implies ( ) = 0. The 

third step is an equality only when \x n \ <C 1; when 02 
is large, e.g. at higher temperatures, or when we use a 
larger free-induction time To to decrease the statistical 
noise level, the variation of x n can be large, and then 
this is not a good approximation. Instead of approxi- 
mating, however, we can compensate the result for the 
sine nonlinearity. Expanding the correlator (Ay m Ay n ), 
we obtain 

(sincc m sinx„) = (sin(u m + v m ) sin(u„ + v„)) = (19) 

= ((sin w m cos ?; m + cosM m smv m )(m — > n)). 

Since sine is an odd function and v n is a zero-mean, Gaus- 
sian i.i.d. variable, (sinu mi „) = 0, and (|19l) becomes 

(sinit m cosu m smu n cosv n ) = (20) 

= (cos v rn ) (cos v n ) (sin u m sin u n ). 

The cosine factors depend on noise in interval 2, i.e., 
above the sampling frequency. This is similar to dephas- 
ing due to quasi-static noise, which acts uniformly on all 
the samples in time (incoherent averaging over a distribu- 
tion of the noise) , and leads to Gaussian decay functions 

(cos cxp(-(T27"o/2). (21) 

For the sine factor, the noise is from interval 1, i.e., it 
is resolved by our sampling, and therefore is not uniform. 
The process is a combination of ensemble-averaged inco- 
herent noise and a frequency-dependent filtering due to 
the (m — n)At time difference in the correlator. Eval- 
uating this factor, we obtain Gaussian damping of a 
hyperbolic-sine function of the correlator, 

(sinu m sinu ra ) = (22) 
= J J du m du n sinu m sinu„ iV(0, o) = 

= exp(-criro ) sinh((u m w n )), 

where the integral is taken over a two-dimensional nor- 
mal distribution N(0,a) with zero mean and correlation 
matrix a = {tr rn , a n , a mn }- (The distribution widths are 
equal, a m = a n , and a mn = (u m u n ) is the correlation 
function.) 
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The correlator (|18|) finally becomes 
(Az m Az n ) = 

= a 2 , expf-cr^o ) expi-a^TQ ) sinh((it m w„)) 



(23) 



Note that no approximation has been made so far (cf. 
Eq. [Toll . If the noise correlation due to l//-type noise 
is small, as in our case where exp(a^Tp) < 10, we can 
neglect the frequency-dependent filtering effect and ap- 
proximate sinh((w m w n )) « (u m u n ). 

Now remains only the determination of the correction 
factors, which we know from the calibration measure- 
ment, exp [(cr 2 + <J2) T o\ = (ao/a(To)) 2 , where we identify 
a(ro)/ao = h(ro) (Eq. [T5|. so that, finally, 



) fa (Az m Az n )/ {a(T )Y 



(24) 



We note that it resembles the signal damping due to de- 
phasing in a Ramsey fringe. The actual numbers used in 
our analysis of the data in Figs. 2 and 4 are presented in 
Tables HfjuU 



C. Data smoothing and reproducibility of the PSD 

The following Figs. [6]-[7] show the reproducibility of our 
results, with sets of data taken on different days. Fig- 
uresEHUshow that our PSD's power laws are independent 
of the choice of smoothing windows. 
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TABLE I: 8e noise (e = 640 MHz). Data in Figs. 2 and 4. 



Temp. (mK) 


to (ns) 


exp(o- 2 To) 


65 


50 


1.5 


120 


50 


1.6 


165 


50 


1.6 


210 


50 


1.8 



TABLE II: SA noise (s = 0). Data in Figs. 2 and 4. 



Temp. (mK) 


ro (ns) 


exp(o- 2 r,?) 


65 


300 


1.2 


120 


300 


1.3 


165 


1,200 


7.4 


180 


1,000 


5.6 



TABLE III: SA noise (e = 0). Data in Fig. 4 (but not in 
Fig. 2). 



Temp. (mK) 


to (ns) 


exp(<T 2 To) 


165 


300 


1.3 


180 


300 


1.6 



FIG. 6: (color online) 5 A noise at 165 mK with two different 
pulse spacings ro, showing reproducibility of the noise-PSD; 
cf. Tables ITlTtilll 
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FIG. 7: (color online) 5 A noise at 180 mK, otherwise like 
Fig. El 
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0.01 0.1 1 10 100 0.01 0.1 1 10 100 0.01 0.1 1 10 100 

Frequency (Hz) Frequency (Hz) Frequency (Hz) 



FIG. 8: (color online) SA noise with different smoothing windows A///. We choose the upper cut-off frequency f c for Fig. 2 
as the lowest frequency for which the phase of the cross-PSD deviates from zero by more than 1 rad. In that figure we use 
the smoothing window A/// = 1/4. The spectrum displays no significant difference depending on A///, and the structure 
can be attributed to insufficient averaging. The phase deviation is, also, due to insufficient averaging, and becomes larger for 
increasing temperature, for a fixed pulse separation To. 




$ 0.01 0.1 1 10 100 0.01 0.1 1 10 100 0.01 0.1 1 10 100 
Frequency (Hz) Frequency (Hz) Frequency (Hz) 



FIG. 9: (color online) Se noise with different smoothing windows A///, cf. Fig. [8] 
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